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ABSTRACT 

We present a photometric study of /-band variability in the young association Cepheus 
OB3b. The study is sensitive to periodic variability on timescales of less than a day, 
to more than 20 days. After rejection of contaminating objects using V, I, R and 
narrowband Ha photometry, we find 475 objects with measured rotation periods, 
which are very likely pre-main-sequence members of the Cep OB3b star forming region. 

We revise the distanc e and age to Cep OB3 b, putting it on the self-consistent 
age and distance ladder of lMavne &: Navlor ( 2008). This yields a distance modulus of 
8.8±0.2 mags, corresponding to a distance of 580±60 pc, and an age of 4-5Myrs. 

The rotation period distribution confirms the general picture of rotational evo- 
lution in young stars, exhibiting both the correlation between accretion (determined 
in this case through narrowband Ha photometry) and rotation expected from disc 
locking, and the dependence of rotation upon mass that is seen in other star forming 
regions. However, this mass dependence is much weaker in our data than found in 
other studies. Comparison to the similarly aged NGC 2362 shows that the low-mass 
stars in Cep OB3b are rotating much more slowly. This points to a possible link be- 
tween star forming environment and rotation properties. Such a link would call into 
question models of stellar angular momentum evolution, which assume that the rota- 
tional period distributions of young clusters and associations can be assembled into 
an evolutionary sequence, thus ignoring environmental effects. 

Key words: accretion, accretion discs, stars:pre-main-sequence planetary systems: 
protoplanetary discs 



1 INTRODUCTION 

There are sound theoretical reasons to expect that accre- 
tion processes help determine the angular momentum of T- 
Tauri stars. Historically, the slow rotation rates of T Tauri 
stars (relative to their break -up velocity) has been explaine d 
by the disc-locking theory IjKoenigl 199ll ; IShu et al. 1994). 
In this theory, magnetic field lines connect the star to 
the disc, enforcing synchronous rotation between the star 
and the material in the disc at some radius, near where 
the magnetic field disrupts the disc. The simplistic the- 
ory has been expanded in a variety of models where angu- 
lar momentum is removed from the star by a combination 
of the disc and an ac c retion-driven wind (e .g. iFendt 20071 ; 
iRomanova et al. 20071 ; iMatt fc Pudritz 20081 ). Whether the 
star is spun up or down by the star-disc interaction depends 



upon the balance between accretion spinning up the star 
and magnetic (or wind) torques slowing it down. Theoret- 
ically, th is issue is unresolve d; some studies find the star 
spins up dBessolaz et al. 2 008). some find it spins down (e.g. 
lLong et al. 20051 ). 

Observationally, the evidence for the influence of ac- 
cretion disks on rotation is mu ch stronger than it was a 
few years ago l|Herbst et al. 20071). Previously, conflicting re- 
sults had arisen (e.g. Herbst et al. 20021 ; IStassun et al. 19991 ; 
iLittlefair et al. 20051 ). These were most likely due to a com- 
bination of small sample sizes, and ambiguous d iagnostics 
of the presence of accretion disks. iRebull et al. f 20051 ) re- 
solved these issues with a large sample of rotation peri- 
ods in the ONC, and accretion disk status defined from 
Spitzer IRAC data. IRebull et al. (2 005) found a clear corre- 
lation between mid-IR excess and rotation, in the sense that 
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stars with mid-IR excess were much more likely to be slow- 
rotators. An ana lysis of the slightly old er NGC 2264 found 
the same result |Cieza fc Baliber 2 007). and confirmed the 
iRebull et al. (20051 ) result, via a refined analysis of the same 
data. There is thus strong observational evidence that the 
star-disc interaction is responsible for extracting angular 
momentum from young stars. Interestingly, a small popu- 
lation of rapidly rotating stars with mid-IR excesses exists; 
it is possible that these stars are being spun-up by accretion, 
hinting that the braking process of young stars is intermit- 
tent. Also, in both the ONC and NGC 2264, there exists a 
significant population of slow rotators with no mid-IR ex- 
cess. These are often interpreted as being recently released 
from disc-l ocking, but there are problems with this interpre- 
tation fsee lBouvier et al. 20071 . for example). 

The firm link established between accretion and rota- 
tion represents significant progress, but there are still open 
questions regarding the rotation of young stars. For exam- 
ple, the low-mass stars are rotating more rapidly than the 
high-mass stars, and appear to spin-up more rapidly as they 
contract toward s the main sequence (jHerbst et al. 2 002; 
llrwin et al. 20071 ). The reason for this is still not known. 
Also, the angular momentum evolution of young stars is de- 
termined by assembling different clusters into an evolution- 
ary sequence, and assuming the period distribution of the 
older clusters can be modelled using the young clusters as a 
starting point. In doing so, the possibility of an environmen- 
tal effect on rotation is ignored. Such an environmental effect 
may be indicated by the data; the young stars in IC348 ro- 
tate much slower than those in the similarly aged NGC 2264 
|Littlefair et al. 20051 ). 

Here we present a photometric study of the young 
association Cepheus OB3b (hereafter Cep OB3b). Cep 
OB3 is a young association covering a region of the 
sky from approximately 22 h 46 m to 23 h 10 m in right as- 
cension and +61° to +64° in declination. The subgroup 
Cep OB3b lies closest to the molecular cloud, and has 
a rich pre-main-sequ ence (PMS) popu lation, confirmed by 
both spectroscopy (|Pozzo et alTljO OS) and X- ray data 
(|Getman et al. 2 006). In section [2] we present the observa- 
tions and data reduction techniques applied. Section [3] de- 
scribes the techniques used to identify periodic variables. 
Section U provides a revision of the age of, and distance to, 
Cep OB3b. In Sections [S] and |§] we present our results and 
analysis of the data, whilst in section [7] we draw our conclu- 
sions. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

RGO /-band CCD images were taken with the Wide Field 
Camera on the Isaac Newton Telescope (INT) on La Palma, 
equipped with 4 thinned EEV 2kx4k CCDs. The total sky 
coverage per exposure is 0.29 sq deg. A single field, with 
the centre of CCD #4 pointing at a = 22 h 55 m 43.3 s , 8 = 
+62 d 40 m 13 3 J2000, was observed, with data being taken on 
every night between 21st September 2004 and 6th October 
2004. This dataset will be referred to as the "short baseline" 
(SB) dataset. Because the SB dataset is not sensitive to 
periods longer than 7-10 days, we augmented it with two 



additional datasets, both of which have longer baselines, but 
less dense sampling. Additional data was taken with the 
wide field camera on 29 nights between 23rd August 2005 
and 1st November 2005. This dataset is referred to as the 
"long baseline 0" (LB0) dataset. 

For the SB dataset the seeing varied between 0.8 
and 5", with a median and standard deviation of l.l"and 
0.9"respectively. Most nights were affected by thin to heavy 
cirrus cloud, although the nights between on 26th September 
2004 and 29th September 2004, and the night of the 5th Oc- 
tober 2004 were photometric. The LB0 dataset was taken in 
seeing conditions between 0.8 and 2.6", with a median and 
standard deviation of 0.9"and 0.8"respectively. Most nights 
were affected by light cirrus only, with less than half a mag- 
nitude effect on transparency. To increase the dynamic range 
of the SB and LB0 datasets, we used exposure times of 5, 
30 and 300 seconds. For the SB dataset these were repeated 
many times throughout a night; for the LB0 dataset one or 
two exposures were taken each night. In total the SB dataset 
contains 477 useable short exposures, 617 useable medium 
exposures and 614 useable long exposures. The LB0 dataset 
contains 41 useable short exposures, 41 useable medium ex- 
posures and 36 useable long exposures. 

The second long-baseline dataset consists of SDSS i! - 
band CCD images obtained with DillCam on the 2-m 
Faulkes Telescope N (FTN). DillCam has a single CCD with 
a pixel scale of 0.278 arcsec per pixel in the default 2x2 
binning mode, giving a 4.7x4.7 arcmin field of view. Four 
fields were observed, chosen to cover the positions of high- 
est stellar density. Data were taken on 15 nights between 
13th October 2005 and 17th December 2005, with between 
2-4 separate observations per night. Exposure times of 30 
and 300 seconds were used. The number of exposures varies 
from field to field, and is different for the short and long 
exposure times, but ranges from 43 to 60 exposures in total. 
This dataset is referred to as the "long baseline 1" (LB1) 
dataset. The LB1 dataset was taken in seeing conditions be- 
tween 1.0 and 3.1", with a median and standard deviation 
of 1.3" and 0.6"respectively. The data are only marginally 
affected by cloud, with transparency variations always less 
than 0.3 mags. 

A plot of the Cep OB3 region and the fields observed is 
shown in figure [1] 

Over the nights of September 12th, 13th and 14th 2003, 
we used the INT to obtain three sets of images of the same 
field as is covered by our /-band data, each set consisting of 
a 300s and a 30s exposure in Ha and a 30s and a 3s expo- 
sure in R. We debiased and flatfielded the images in a similar 
manner to the /-band data, but then instead of searching 
for stars in a combined image, we carried out optimal pho- 
tometry at the posi tions of all the objects in the catalogue of 
iMavne et al. f2007l ). See Section \2. 21 for a complete descrip- 
tion of how the photometry was carried out. We then profile 
corrected an d combined the mag nitudes from each image, as 
described in lMavne et al. f2007l ) , to obtain R and R — Ha 
magnitudes, with an arbitrary zero p oint. We then combi ned 
this with the V and I photometry of IMavne et al. f2007l ) to 
create Table 1 (the full version of which is available elec- 
tronically). 
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Table 1. Photometry for the INT WFC field observed is available in this format at 
http://www.astro.ex.ac.uk/people/timn/Catalogues/tables.html This table is a sample only, to indicate the format of the online table. 
The flag column indicates data quality, and is explained in full at http://www.astro.ex.ac.uk/people/timn/Catalogues/format.html 



Field/CCD (field.ccd) ID RA (J2000) Dec (J2000) xpos (ccd) ypos (ccd) MAG uncertainty flag COL uncertainty flag 



1.02 410 22 52 47.929 +62 30 57.08 613.729 2634.474 11.723 0.010 OO 0.763 0.014 OO 

1.04 10 22 57 13.885 +62 41 35.55 774.439 119.826 11.596 0.010 OO 0.589 0.013 OO 




Figure 1. INT WFC field-of-view overlayed on DSS2 red image 
of Cep OB3. Also shown are the four Faulkes telescope point- 
ings obtained. Circles denote the locations of periodic association 
members (see section \3\ for details). 



2.2 Image Processing and Optimal Photometry 

Image processing and data reduction was performed in a 
consistent manner for all datasets. The individual frames 
were bias subtracted using a median stack of several bias 
frames. The frames were then flat fielded using twilight sky 
flats taken on one of the photometric nights. Only sky flats 
with peak counts of less than 30,000 were used, to avoid 
non-linearity effects. A bad pixel mask was constructed for 
all frames by flagging all pixels which deviated by more than 
10(J from the median of a ratio of two flat fields. In addition, 
a smoothed version of a long exposure was subtracted from 
the long exposure, and all pixels with a value lower than 
1000 were masked. These procedures accounted for most bad 
pixels, but left a number of bad columns unmasked. These 
columns were identified by hand and added to the bad pixel 
mask. 

The processed CCD frames were then analysed using 
optimal photometry, as im plemented by the c luster soft- 
ware described in det ail bylNavlor et al. (20021 ). with modi- 
fications described bv lLittlefair et aL (20051 ). Advantages of 
this approach over classical aperture photometry include 
better signal-to-noise ratios, and robustly determined un- 
certainties for each observation. Photometric measurements 
which are deemed of poor quality are flagged, rather than 
rej ected. The reduction of our datasets closely follows that 
of iLittlefair et al. (20 051. Profile corrections (analogous to 
aperture corrections in classical aperture photometry) were 



allowed to vary as a 2nd order polynomial function of posi- 
tion upon the CCD. 

After profile correction, the photometric measurements 
were adjusted for any difference in the airmass and trans- 
parency for each frame, by determining a relative trans- 
parency correction from the bright stars. Before this pro- 
cess we added an additional, magnitude-independent er- 
ror of 0.01 mags to the results of each frame, in order to 
yield a plot of %v versus signal-to-noise ratio that was flat 
and had a modal value of approximately 1. In the absence 
of colour information, transparency correction may intro- 
duce errors into our relative photometry because of the 
colour dependence of extinction. Following the analysis of 
ILittlefair et al. (20051 ). we find that any errors introduced 
should be smaller than 3mmags. Hence we are confident 
this effect is negligible for the purpose of our analysis. An 
astrometric solution was achieved through comparison with 
a 2MASS catalogue of the same region. A 6-coefficient so- 
lution to yielded a rms discrepancy in positions which was 
always less than 0.14 arcseconds. 



2.3 Transformation to a standard system 

Our aim for this dataset was to obtain high standard rel- 
ative photometry of Cep OB3b. This meant that we did 
not measure colour information for each star at each epoch, 
and hence we are unable to tie our photometry to a stan- 
d ard system directly. Instead, we use the BVI photometry 
of iMavne et al. (20071 ) throughout this paper. 



2.4 Final Dataset 

The final data set consists of lightcurves for 42962 stars. 
Rather than combine data with very different sizes of error 
bars, three lightcurves were produced for each star, resulting 
from the 300, 30 and 5 second exposures (where available). 
This was done for each of the SB, LBO and LB1 datasets, re- 
sulting in a maximum of 8 lightcurves for each star. Figure[2] 
shows the RMS variability in magnitudes for our stars, plot- 
ted as a function of magnitude. This plot shows the high 
internal accuracy reached in our dataset. In both the SB 
and LBO datasets, obtained with the INT WFC, the 300s 
and 30s datasets have an internal accuracy of better than 1 
per cent, whilst the 5s second dataset has an internal accu- 
racy of 2 per cent. The gradual decrease in internal accuracy 
with decreasing exposure time is most likely due to a corre- 
sponding decrease in the number of stars available to per- 
form the transparency correction. The LB1 dataset, taken 
on the FTN, reaches an accuracy of 1 per cent in both the 
30 and 300s exposures. 
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Figure 2. Scatter in the photometry of stars as a function of 
brightness in the SB (top), LB0 (middle) and LB1 (bottom) 
datasets. The 300s exposures are plotted in red (dark grey), the 
30s exposures in green (light grey) and the 5s exposures in blue 
(black). The dotted line marks a photometric precision of 1% 

3 DETECTION OF PERIODIC VARIABLES 

To search for periodic variables within our datasets, Lomb- 
Scargle periodograms were calculated independently for 
each of the SB, LB0 and LB1 datasets. Frequencies corre- 
sponding to periods from 0.1-12 days were searched for the 
SB dataset, and 4-25 days for the LB0 and LB1 datasets. 
In each case the lower value is set by the "mean nyquist" 
frequency of the lightcurve. 



Searching for periodicities in photometry of young 
stars is a challenging task. The unevenly spaced data 
sampling means that analytical estimates of the false 
alarm probability (FAP) ( e.g. Irlorne fc Baliunas 19861 ; 
iLinnell Nemec fc Nemec 19851 ) are rendered invalid. Worse 
still, many young stars show irregular variability as well as, 
or instead of, periodic photometric modulations from sur- 
face spots. The interaction of this intrinsic variability with 
the often patchy sampling window can quite easily introduce 
peaks into a periodogram which appear significant based on 
a simple FAP cut. 

The inadequacy of analytical estimates of FAP is of- 
ten sidestepped by estimating the FAP via a Monte-Carlo 
simulation; the FAP is set to be the fraction of simulated 
lightcurves where the peak power exceeds the observed peak 
power. The simulated lightcurves can either have pure Gaus- 
sian noise added, or some correlated noise model can be as- 
sumed. Eventually, one ends up with a cut-off level for the 
peak power, above which a period is judged to be significant. 

The goal of these Monte-Carlo simulations is to estimate 
the heights of spurious peaks which may arise in the peri- 
odogram due to a combination of imperfect sampling, and 
variability in the lightcurve. However, they can only realis- 
tically account for variability in the lightcurve due to pho- 
tometric error. For a large number of young stars, intrinsic 
variability is larger than the photometric errors. In addition, 
colour-dependent effects mean that transparency variations 
can be imperfectly corrected, introducing a night-to-night 
variability that Monte-Carlo simulations don't attempt to 
account for. Selecting periodic variables solely on the ba- 
sis of Monte-Carlo simulations to asses FAP can thus lead 
to significant contamination of a dataset by spurious peri- 
ods. The large number of objects with periods very close 
to Id in the NGC 2264 dataset of lLamm et al. (20051 ) and 
lLamm et al. (20*53 ) may be an example of this. 

What we are aiming for is a sample of genuine pe- 
riodic variables, which is relatively free from contamina- 
tion by spurious periods. Since we have a large sample of 
lightcurves to search, one approach is to use the lightcurves 
themselves as a means of estimating the height of spuri- 
ous peaks introduced into the periodogram by all factors 
affecting the dataset, from imperfect photometry to intrin- 
sic variabihty ; _We_ca^ analytical FAP, as estimated 
bv lHorne fc Baliunas (19861 ). FAPh, for the strongest peak 
in all our lightcurves. Assuming that most stars in our field 
of view are not periodic, the distribution of FAP;, is de- 
termined by the interplay of stellar variability, finite sig- 
nal/noise and systematic effects in our photometry with the 
sampling windows of our lightcurves. Candidate periodic 
variables can be selected as outliers in the FAP;, distribu- 
tion. The value of FAPj, to make a cut at is selected from 
a histogram of the FAPj, values from our lightcurves. The 
value of FAPh is chosen by-eye from the histograms shown 
in figure [3] After applying the additional selection criteria 
described below, the level of contamination is assessed by vi- 
sual inspection of randomly selected folded lightcurves, and 
the chosen FAPh value is adjusted to keep the contamina- 
tion at an acceptable level. During the visual examination 
of folded lightcurves, most rejected candidates were those 
whose periodic nature rested on only a few, discrepant points 
(such as might occur from two flares repeated days apart, 
for example). Lightcurves were inspected by a single author 
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with poor rotational phase coverage. Repeating patterns in 
bright, irregular variables can add significant power to a 
Lomb-Scargle periodogram, and stars with poor phase cov- 
erage must necessarily be considered dubious periodic can- 
didates. To remove these interlopers we applied two further 
cuts: one on the \v with respect to a sinusoidal fit at the sug- 
gested period, and another on a stat istic designed to max - 
imise phase coverage, the S-statistic l|Saunders et al. 20061 ). 
The S-statistic is the sum of the squares of the phase dif- 
ferences between adjacent data points, after the data has 
been sorted according to phase. It is normalised by division 
by the number of datapoints, such that uniform phase cov- 
erage gives an S- statistic of unity, and the S-statistic rises 
with increasingly irregular phase coverage. The S-statistic 
is especially good at removing spurious periods arising from 
the 1 day -1 natural frequency that is introduced into the 
data the rotation of the Earth, and observing at a single 
site. The S-statistic for our datasets as a function of period 
is shown in figure [4] This figure shows the improvements in 
phase coverage offered by the long baseline datasets; partic- 
ularly for periods around multiples of a dajQ. We calculate 
the S-statistic for all the periodic candidates in our sample, 
after converting to phase using the best fitting period as 
deduced from the Lomb-Scargle periodogram. We rejected 
all stars with an S-statistic greater than 30 (SB), 5 (LB0), 
and 7 (LB1), or with periods in the ranges 0.9-1.05 d, 0.495- 
0.505 d and 0.325-0.335 d. The X l cut was established via a 
by-eye inspection of 200 randomly selected candidate peri- 
odic variables. The value chosen rejected erratic variables, 
whilst preserving the majority of stars which were clearly 
periodic, but non-sinusoidal. We are thus confident that the 
Xt cut is not rejecting periodic variables with a high sig- 
nal/noise ratio whose lightcurves are not purely sinusoidal. 
All stars with \ 2 V > 1-3 were rejected. 

The selection process outlined above is far from per- 
fect, but there is no ideal solution yet found to the problem 
of searching for periods in an imperfect dataset of objects 
which exhibit significant intrinsic variability. The selection 
criteria adopted above are conservative; the aim was not to 
discover all rotational periods, but to produce a catalogue of 
periodic stars which was relatively free from contamination 
by spurious periods. To check our success in this aim we 
undertook a by-eye examination of 200 randomly selected 
folded lightcurves. This examination suggests that our final 
sample of periodic stars has a contamination by spurious 
periods of ~3%. 

Once each dataset was searched for periods, we com- 
bined the results to produce our final catalogue of periodic 
variables. As each star could have a detected period in each 
of the datasets, we selected a period from the dataset which 



Figure 3. The distribution of FAP h , for the SB (top), LB0 (mid- 
dle) and LB1 (bottom) datasets. 



(Littlefair), and the process repeated until the results were 
acceptable. We settled on a FAP h cutoff of 10 -35 for the SB 
dataset, and a FAPh cutoff of 0.005 for the LB0 and LB1 
datasets. 

A FAP cut taken alone will not reject some objects 
which should not be included in a sample of periodic vari- 
ables. Two examples are very irregular variables and stars 



1 The long period datasets show spikes in the Saunders statistic 
at periods between 11 and 15 days. This means the phase cover- 
age of these datasets degrades when the data is folded on these 
periods. In figure [5] we show a graphical representation of the 
temporal coverage of each of our datasets. Looking at this figure, 
we can understand that the spikes in Saunders statistic at 11-15 
days arise because the long period datasets have a tendency to- 
wards gaps in coverage lasting 5-8 days. We are not able to offer 
a definitive explanation of the reason for these gaps in coverage, 
but we can speculate that it is linked to the typical duration of 
bad weather at our observing sites. 



6 



Long Baseline 1 



I 

i 0.5 ■ 



Integer part of date (days) 



Long Baseline 




1 



3D 35 




3D 35 



Short BoseLmf- 



Figure 6. The fraction of correctly recovered periods, as a func- 
tion of period, for a representative lightcurve from each dataset. 
The red histogram shows the sensitivity to sinusoidal periods with 
signal/noise ratios of 5, the grey histogram shows the sensitivity 
to sinusoidal periods with signal/noise ratios of 2. 
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Figure 5. Time sampling of our datasets. The x-axis shows the 
integer part of the time elapsed since the first observation. The 
y-axis shows the fractional part of the time elapsed since the first 
observation. 



was most likely to be accurate. For periods longer than 7 
days we adopted a period from the dataset with the longest 
baseline, whilst for shorter periods we adopted a period from 
the dataset with the densest time sampling. Where a dataset 
has lightcurves with differing exposure times we always used 
the deepest lightcurve in which a period was detected. In 
the vast majority of cases, periods from different datasets 
were consistent. Only 8 objects showed significant disagree- 
ment between datasets; these were long-period systems for 
which the SB dataset lacked sufficient baseline. The final 
periodic dataset consists of 709 periodic variables; 578 from 
the SB dataset, 88 from the LBO dataset and 43 from the 
LB1 datas et. 704 of our period ic variables have BVI photom- 
etry from lMavne et al. (20 07). These 704 objects comprise 
our final catalogue of periodic variables, and are presented 
in Table 2, which is available in the online version of this 
paper. 



3.1 Period completeness and selection effects 

We investigated our ability to recover a range of periods 
by introducing artificial sinusoidal signals into our data 
and checking to see if the correct period was recovered. 
Lightcurves from our datasets were selected at random, nor- 
malised and had an artificial sinusoidal signal added with a 
signal-noise of 2. For our dataset, this is equivalent to adding 
a 0.02 mag signal to an object with 1=18 (the faintest mag- 
nitude for which report periods, and corresponding roughly 
to a mass of 0.2 Mq). Our simulation is thus directly com- 
parable to the simulation bv llrwin et al. (2008) for their pe- 
riodic dataset in NGC 2362. We also performed a simulation 
in which sinusoidal signals with a signal-noise ratio of 5 were 
added to a normalised lightcurve. 

The period selection process described above was then 
performed on our artificial lightcurves. A period was deemed 
correctly recovered if the if the recovered period was within 
10 per cent of the injected period. Figure [6] shows the re- 
sults. The SB dataset finds essentially all periods below 7 
days, but is incomplete above that. This is expected from 
the relatively short baseline of this dataset; indeed the rea- 
son so few periods are recovered above 10 days is that these 
periods fail our test based on the Saunders statistic. The 
LBO dataset improves upon this substantially, but is still 
not complete above 10 days. Because of the relatively sparse 
sampling of this dataset, completeness for the LBO dataset 
is a strong function of the injected signal size. For weak sig- 
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Table 2. Data on the detected periodic variables is available in this format at http: / /www. astro. ex. ac.uk/people/timn/Catalogues/tables. html 
This table is a sample only, to indicate the format of the online table. The flag column indicates data quality, and is explained in full at 
http: / /www. astro. ex. ac.uk/peoplc/timn/Catalogues/format. html 
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Figure 4. The Saunders statistic, as calculated from a representative lightcurve from each dataset, in which more than 60 per cent of 
the data points in the lightcurve were unfiagged. 



rials the completeness never rises much above 50 per cent, 
whereas this dataset recovers most stronger periodic signals 
below 10 days. The LB1 dataset is largely complete at all 
periods, with the exception of a small drop in completeness 
between 12-14 days, however, the small spatial coverage of 
this dataset means that it will not be available for many 
stars. 

This simulation shows that all datasets are sensitive to 
short periods (below 7 days). In practise, the SB dataset will 
be much better at recovering period shorter than 7 days. 
This is because the simulation described above does not 
take account of the stochastic variability often seen in young 
stars. This variability can be considered red noise, and has 
characteristic timesca les of days to weeks. As discussed by 
iLittlefair et al. (2005h . dense sampling of lightcurves makes 
it much easier to recover periods in the face of this red noise. 
Therefore, we prefer periods from the SB dataset for periods 
less than 7 days. 

In figure [7] we show the X-ray counts from 



iGetman et al. (20061 ) for those objects for which we found 
periods, and for those objects in which no periods were 
found. Objects with detections in the X-ray catalogue of 
Gct man et al. (2006) were matched with our optical cata- 
logue using a 3 arcsecond matching radius. A cut in colour- 
magnitude space (see section |3~2")| was applied to our optical 
catalogue to reduce the contamination from foreground ob- 
jects. The matches were then divided into two groups; those 
in which we detected periods, and those in which no period 
was found. We see a suggestion that the stars in which we 
found periods are significantly more X-ray luminous than 
those stars for which a period was not found. This result 
is not formally significant; A 1-D K-S test gives a 13 per 
cent chance that they were drawn from the same parent dis- 
tribution. Thus it is possible, though not certain, that our 
catalogue of periodic objects is biassed against X-ray faint 
sources. This bias is presumably introduced because strong 
starspot activity is associated with increased X-ray emission. 
Whether this introduces bias into our period distribution 
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Figure 7. The histogram of X-Ray counts for objects with known 
periods (bottom panel) and for objects with no known periods 
(top panel). The X-Ray flux of objects for which we found periods 
is systematically higher than for objects where no period was 
found. 



remains to be seen. There is no evidence for a correlation 
between X-ray counts and rotational period within the stars 
for which periods were found. 

To summarise; our survey is sensitive even to weak pe- 
riodic signals below 7 days and down to masses of 0.2 M©. 
For longer periods, completeness is somewhat lower, and we 
are, generally speaking, only sensitive to stronger periodic 
signals (S/N ^ 5). Thus, our periodic dataset is biassed to 
some degree towards shorter rotational periods. 

3.2 Rejection of background and foreground 
objects 

There will also be some contamination of our sample by 
non-PMS stars. Whilst variability is a very good indica- 
tor of youth, periodic variables do exist amongst the field 
population, and these should be removed from our sample. 
The most practical way of removing foreground and back- 
ground objects is via cuts in colour-magnitude space. PMS 
objects are larger and correspondingly brighter than main 
sequence stars of a given colour and distance. A simple cut 
in colour-magnitude space can remove a large amount of 
the contamination by main sequence and other foreground 
objects. We decided on an appropriate cut by-eye, which 
is shown in figure [8] Of the 704 periodic stars, 553 have 
unflagged colours placing them above and to the right of 
the cut in figure fS] Our dataset is still not free of possible 
contamination. lLamm et al. (20*041 ) points out that signif- 
icant contamination from variable background giants (e.g. 
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Figure 8. The colour magnitude diagram of periodic variables 
in Cep OB3b. An initial colour cut, used to select PMS objects 
is plotted with a dashed line. 



RR Lyrae stars) can lie within the PMS region as defined 
by the cut a bove, and a further cut is required to remove 
these objects. "Lamm et al. (20041 ) suggest that the difference 
in (R — Ha) colours between giants and main-sequence or 
pre-main-sequence stars allows for rejection based upon this 
colour index. In addition, the extra reddening experience by 
background giants gives them larger (V — I) colours than 
main-s equence or pre-ma in-sequence stars. Therefore, fol- 
lowing |Lammelal~(200l), we use the (R — Ha) vs. (V — I) 
colour-colour diagram, shown in figure [5] to reject back- 
ground giants. The locus of PMS stars was defined as the 
median (R ~ Ha) colour as a function of (V — I) for all pe- 
riodic stars which survived our initial cut in V vs. (V — I). 
A quadratic fit to the median defines the PMS locus as 



(R - Ha) loc 



-0.024(1/ - If + 0.179(1/ - /) - 3.126(1) 



To reject background giants we define a lower cut-off, shown 
in figure [9] as a dashed line, by (R — Ha) — 1.655, where S 
is a fit to the standard deviation in (R — Ha) as a function 
of (V — I), given by 



S = exp [0.66(1/ - /) - 4.23] . 



(2) 



Of the 553 stars which passed our V vs. V — I cut, 
481 had unflagged colours in both (R — Ha) and (V — I), 
6 of which were rejected as background giants. Our final 
dataset thus consists of 475 periodic variables which are very 
likely pre-main-sequence members of the Cepheus OB3b 
star-forming region. 



4 AGE AND DISTANCE OF CEP OB3B 



In iMavne et al. (20071 ) we derived ages for a group of 
well studied star-forming regions, including Cep OB3b, 
based on the fact that the luminosity of the pre-main- 
sequence declines with age. Thus the derived age clearly 
depends on the distance assumed for each group, which 
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Figure 9. The (R — Ha) vs. (V — /) colour-colour diagram for 
Cep OB3b. The red (light grey) points represent all the stars 
for which we have photometry. The black dots are those periodic 
stars which survived the cut in the V vs. V — I cmd. The pre- 
main-sequence locus is plotted as a solid line, and the locus used 
to reject background giants is plotted as a dotted line. 
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Figure 10. The members of Cep OB3a (circles with error bars) 
and Cep OB3b (st ars with error bars). The photometry (from 
iBlaauw et al. 19591) . and the colours have then been de-reddened 
and the magnitude corrected for extinction by deriving colour 
excesses in the U — B vs B — V plane. The magnitudes have the 
been converted into absolute magnitudes using the best-fitting 
distance modulus of 8.8. The curves are Geneva-Bessell isochrones 
for (from faintest to brightest) 3Myr, lOMyr and a 3Myr equal- 
mass binary sequence. 



in iMavne et al. (20071 ) we took from the literature. Real- 
ising these literature distances were the major source of un- 
certainty in the ages, IMavne fc Navlor (20081 ) derived con- 
sistent distances with robust uncertainties by fitting the 
main-sequence stars in each region using the r 2 method of 
iNavl or fc Jeffries (20061 ). We have applie d the same methods 
to the 7 Vel a ssociation and Vel OB2 (jjeffries et al. 20*091 ) 
and NCC2169 dJeffries et al. 20071). which along with the re- 
gions presented in Mavne fc Navlor (20081 ) represent a self- 
consistent set of ages and distances for a significant group 
of young clusters and star-forming regions. 

Clearly we must attempt to fit Cep OB3b into 
this system, but a lthough we attempted this in 
IMavne fc Navlor (20081 ). we failed to do so because we 
had only a few stars which could be used for main-sequence 
fitting, and they gave an answer sharply at variance with 
the literature values. In what follows we use the stars from 
Cep OB3a as well as Cep OB3b, and still obtain an answer 
which is different from that determined by others, but we 
now un derstand that the diffe rence is due to the earlier work 
using a Sc hmidt-Kaler (19821 ) main-sequence (as opposed to 
his zero-age main-sequence), which is significantly brighter 
than more modern work. 

To determine an age in the same way as 
IMavne et al. (20071 ) we must place the pre-main-sequence 
of Cep OB3b in a Vb vs ( V — I) a CMD, and compare it with 
similar sequences for other star-forming regions. If we carry 
out such a comparison by plotting the members of each 
group in a CMD the number and spread of the data points 
m ake the plot impossi ble to interpret. Instead, as described 
m IMavne et ai~(2 007). we use splines fitted through the 
members. To place Cep OB3b in such a diagram we 
need the apparent distance modulus (to correct V), and 
E(V — I) to correct V — I. In what follows we derive these 



by dereddening an d then fitting the main sequence, in an 
identical fashion to lMavne fc Navlor (2 008). 



4.1 MS star Individual reddenings and extinctions 

iPozzo (2001) provide d a list of members of both sub-groups 
with photometry in (|Blaauw et al. 19591 . hereafter BHJ).0 
We dereddened each of these in U — B vs B — V and used 
the implied extinction in V to create the CMD shown in 
Figure 1101 To do this we have used the extinction vec- 
tors of iBessell et al. (1998a), and a Geneva-Bessell main- 
sequence (see | Leieune fc Schaerer 200ll ; lBessell et al. 1998al ; 
IMavne fc Navlor 2008T T The mean extinction of the sample 
is E(B -V) = 0.79 with an RMS of 0.16 mags. 



4.2 True distance modulus 

For our distance modulus determination we must ensure we 
fit only those st ars which are on the m ain sequence (see 
the discussion in IMavne fc Navlor 20081 ). The bright limit 
to the main sequence is defined by the turn-off, which be- 
comes fainter with age. In Figure [10] we show the 3 and 
lOMyr isochrones, which make it clear that the turn-off is 
not a sharp transition. We chose these ages as 3Myr is the 
youngest age available in the Geneva isochrones, whilst it 
is unlikely that any members of Cep OB3 are older than 10 
Myr. We therefore fit only those stars lying below V = —1.5, 
since brighter than this the two sequences separate by more 



2 The star numbers in BHJ for these are as follows. Sub-group a: 
14, 18, 19, 23, 37, 44, 46, 47, 50, 54, 56, 59, 66, 68, 69, 70, 75, 76, 
77. Sub-group b: 2, 8 , 10, 11, 15, 16, 17, 20, 22, 24, 25, 26, 31, 
33, 39, 40, 41 
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Intrinsic B-V 

Figure 11. The fit to the main-sequence data. The data are 
shown as circled error bars , and the colour scale is the best fitting 
model (p in Equation 3 of lNavlor &: Jeffries 200(|). 

than 0.009 mags in B — V, correspon ding to half the uncer - 
tainty for the photometry quoted bv lBlaauw et al. f 19591 ). 
The faint end of the main sequence, within the CMD, is ap- 
proximately defined by the radiative-convective (RC) gap, 
which is only defi ned (and observabl e) for clusters younger 
than 20 Myrs old IjMavne et al. 20071 ). Although some stars 
blueward of the RC gap will not have reached the main- 
sequence, their position in CMD space is very close to 
the main-sequence loc us at these ages (see discussion in 
iMavne fc Navlor 2008h . Even at the age of the Orion Neb- 
ula Cluster (2Myr) the stars bluer than B — V — are on 
the main sequence, so the entire remaining sample can be 
fitted. 

We can now fit a main-sequence isochrone to the re- 
maining d ata in Figure | 10| using th e r 2 fitting procedure de- 
sc ribed inlMavne fc Navlor (20 08). with the improvements 
m iNavlor (2009^ Using the Geneva-Bessell isochrones we 
obtain the fit shown in Figure 1111 and a true distance mod- 
ulus of 8.8 ± 0.2mags. 

4.3 Apparent U-band distance modulus 

Since our aim is to correct the PMS data of 
IMavne et al. (20071 ) into the Vo vs (V - I) CMD, we 
need the apparent distance modulus for the field from 
which those data were taken. This is the same INT field 
as observed here, within which there are twelve BHJ stars, 
which give a mean E(B — V) = 0.83 ±0.11 where the 
quoted uncertainty is the RMS about the mean. We can 
then calculate Ay as 3.26 x 0.83 = 2.71, and the apparent 
F-band distance modulus as 8.8 + 2.7 = 11.5 ± 0.2. This 



is independent of any assumptions about R = ' 
since changing R changes the extinction-free magnitudes of 
Figure [10] and hence the true distance modulus by an equal 
an opposite amount to the change in U-band extinction. 



4.4 Reddening in V — I 

Finally we need E(V — /) to correct V — I. The usual 
way of proceeding at this point would be to use the ratio 
E(V - I)/E(B - V) to calculate E{V - I) from E(B - V). 
But in this case there is a subtle issue which must be 
addressed. For a given column of interstellar material be- 
tween the observer and a star the extinction (e.g. A v ) and 
reddenings (e.g. E(B — V)) depend on the colour of the 
star being observed. Therefore to specify the column den- 
sity precisely, one must specify the colour of the star for 
which the measurement was made. For B-stars let us re- 
fer to E(B -V) B , and for the red PMS stars E(B - V)r. 
To construct the normal colour dependent reddening vec- 
tors one constructs a continuous sequence of ratios from say 
E(V - I)b/E(B - V) B to E(V - I)r/E(B - V) R . Our 
case is somewhat different, we have measured E(B — V)b 
but wish to know E(V — I)r, and therefore need the ra- 
tio E(V - I) R /E(B - V)b- This can be obtained folding 
a blue stellar atmosphere through the appropriate filter re- 
sponses, and then applying a given column density to the 
flux, to obtain E(B — V)b- Using the same column density 
one can also obtain E(V — I)r. Whilst iBessell et al. (1998al ) 
find that E(V - I)r/E(B - V) R = 1.44 for (V - 7) = 2, 
we fin d that with identical atmospheres and the extinction 
law of lCardelli et al. (19891 ) E(V - I) R /E{B -V)b = 1.26 
for red stars of (V — I)q of approximately two. 

In principle Ay should show a similar effect, with red 
stars being less affected by a given column density than blue 
ones, and so, in principle we have overestimated the U-band 
extinction for the red stars. However, we have calculated 
this effect, and find it is less than 0.03 mags, and so ignore 
it. Thus our final estimate of E(V — I) for appropriate stars 
of PMS colours is 1.26x0.83=1.05. 



4.5 Age estimate 

Now we have distance moduli we can obtain a relative age for 
Cep OB3b by comparing the absolute magnitude of its pre- 
main-sequence with those of other young groups. Within any 
young (<5Myr) OB association or cluster, there is a large 
scatter in the luminosity of the pre-main-sequence mem- 
bers at a ny given colour (for a d iscussion of this see, for 
example, iBurningham et al. 20051 . and references therein). 
To make the comparison between groups easier, therefore, 
we fit a curve through the members of each group, and 
then compare the positions of these fits in an absolute- 
magnitude intrinsic -colour diagram. The fits are taken from 
IMavne et al. (20071 ) and were constructed by binning the 
data by magnitude, finding the median magnitude and 
colour within each bin, and then fitting a spline through the 
resulting points. The corrections to abs olute magnitude and 
in trinsic colour were taken from either IMavne et al. (20071 ) 
or lMavne fc Navlor (20081 ). except for Cep OB3b for which 
we used the values derived above. The result is shown in Fig- 
ure [T2j which shows the position of Cep OB3b very closely 
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Figure 12. A comparison of the position of the PMS in a CMD 
of Cep OB3b with other young groups. The solid curve is for 
Cep OB3b. The dashed curves are the sequences for (from left 
to right) the Orion Nebula Cluster (2Myr), NGC2362 (lower, 4- 
5Myr), IC348 (upper 4-5Myr) and h and \ Per (13Myr). 



matches that of IC348 and NGC2362, both which were as- 
signed ages of 4-5Myr, and so we will use an age of 4.5Myr 
for Cep OB3b throughout this paper 



4.6 Comparison with literature values 

In iMavne et al. (20071 ) we adopted an age of about 3Myr, 
and a distance modulus of 9.65 ±0.2 mags. Our revised true 
distance modulus of 8.8 m ags is much lower, ye t our new age 
is only slightly older. The lMavne et al. (20071) distance wa s 
essentially the one derived by iMoreno-Corral et al. (19931 ). 
which relies on the ISchmidt- Kaler (19821 ) main sequence for 
stars around the solar neighbourhood. This is about 0.8 
mags brighter in Mv than the Geneva-Bessell isochro ne for 
3Myr, and indeed the zero-age ISchmidt-Kaler (19821 ) main 
sequence. This effect alone would push Cep OB3b to an age 
of around lOMyr, but is partially compensated for by our 
value of E(V - I)/E{B - V) . The difference betw een the 
value we use, and say that of iBessell et al. (1998al ). results 
in shift of 0.15mags in V — I, decreasing the age from lOMyr 
to the 4.5Myr adopted here. 

Given the significant change in the distance compared 
with literature values, and the low value of the V — I colour 
excess we have used, we have checked that the resulting age 
is consistent with other available age indicators. These are 
as follows. 

(i) Cep OB3b has a large R-C gap (see IMavne et al. 20071 ) 
suggesting an age of around 3Myr. 

(ii) The PMS shows a large "age spread" in the CMD, which 
IMavne et al. (20071 ) show is characteristic of groups younger 
than about 5Myr. 

(iii) The turn-off age derived by iJordi et al. (1996) is 
5.5Myr, and although revised distance modulus would make 
this slightly younger, the turn-off and PMS ages are clearly 
not grossly discrepant. 

(iv) There is still molecular material around the group (it is 
obvious in the optical images and the extinction is patchy), 



which is again normally associated with young groups. Thus 
we conclude that all current age estimates are in concor- 
dance. 



5 PERIOD DISTRIBUTION 

5.1 Mass dependence of rotation 

The dependence of ro tation period upon mass was first 
shown for the ONC b vlHerbst et al. (20021). a nd later con- 
firmed in NGC 2264 l|Lamm et al. 200411 20051 ) and IC 348 
jLittlefair et al. 20051 ). In general, lower mass stars rotate 
significantly faster than higher mass stars, both in the 
PMS phase and in you ng clusters (e.g. iHerbst et al. 20071 ; 
IScholz et al. 2005l f 20091 ) . The reasons for this are not yet 
clear, though hints are emerging that it is related to a 
change in magneti c structure with mass l|Scholz et al. 20091 ; 
iDonati et al. 20091 ). 

We calculated masses for our stars by comparing ex- 
tinction and distance corrected 7- band m agnitudes with the 
4.5Myr isochrones of lBaraffe et al. (19981 ). For this we need 
the apparent distance modulus in 7, which we calculate from 
the true distance modulus (Section 14. 2 1 and Aj. We do this 
using the ratio Ai/E(B — V)b, calculated for red stars in 
the same way as descri bed in Sectio n 14.41 We find this ratio 
is 1.94 (c.f. 1.89 from IBessell et al. 199 8b). giving an ap- 
parent 7-band distance modulus for Cep OB3b of 10.4. We 
also calculated masses for the periodic stars in NGC 2264 
and NGC 2362. We adopted ages, true d istance moduli and 
E{B-V) B from lMavne fc Navlor (20081 ). Ai is calculated as 
outlined above to obtain the apparent distance modulus in 
7, and the masses a re derived by comparison with the same 
iBaraffe et al. (1 9981 isochrones at the appropriate age. For 
the ONC a similar method was use d, except that indiv idual 
de-reddened 7 band magnitudes of iHillenbrand (1 997) were 
used. In this way, we have ensured that all clusters have 
masses and ages which are determined in a self-consistent 
manner. The results are shown in Figure 1131 which shows 
the rotation periods of young stars in Cep OB3b, NGC 2362 
NGC 2264 and the ONC, plotted as a function of mass. 

We can see from Figure Q3] that in many young clus- 
ters there is a strong dependence of rotation with mass, in 
the sense that the stars with masses below 0.4 Mq show a 
dearth of slow rotators. This can readily be seen in the pe- 
riod distributions of NGC 2264 and NGC 2362, for example. 
To study this effect in Cep OB3b, we divide our data into 
high and low mass samples. The high mass sample contains 
stars with masses between 1.0 and 0.4 Mq, whilst the low 
mass sample contains stars with masses between 0.4 and 0.2 
M 0. The rotation distribution in Cep OB3b supports the 
general picture seen in other clusters; the low mass stars 
are rotating, on average, more rapidly than the higher mass 
stars. The median period of the low-mass stars is 3.7 days, 
against 4.4 days for the high mass stars. Unlike in other 
clusters, however, this result is barely statistically signifi- 
cant. A ID K-S test gives a 6 per cent chance that the low 
and high mass stars were drawn from the same period distri- 
bution. Under the assumption that the period distributions 
in high and low-mass stars differ only in their median value, 
a Mann- Whitney test gives a 3 per-cent chance that the me- 
dian period is the same for the high and low-mass stars in 
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It is apparent that the low-mass stars in Cep OB3b are rotat- 
ing more slowly than those in other clusters. A comparison 
with NGC 2362 is most interesting, as it is at a similar age 
to Cep OB3b. The median period for low-mass stars in Cep 
OB3b is 3.7 days, compared to 2.3 days for NGC 2362. A 
ID K-S test gives a 0.1 per cent chance that the periods 
of low mass stars in the two clusters were drawn from the 
same period distribution. Amongst the high mass stars the 
picture is different, with the high mass stars in Cep OB3b 
rotating more rapidly than their counterparts in NGC 2362. 
The median period for the high-mass stars in Cep OB3b is 
4.4 days, compared to 4.9 days for NGC 2362. This differ- 
ence is only weakly significant though; a ID K-S test gives a 
6 per cent chance that the periods of high mass stars in the 
two clusters were drawn from the same period distribution. 
It is not clear why the two clusters should be so different. 
Our simulations (see section I3.1|l suggest that we should 
have found almost all periods below 7 days with amplitudes 
greater than 0.02 mags in stars down to 0.2MQ. Similar cal- 
culations by llrwin et al. (20081 ) suggest that the survey in 
NGC 2362 was sensitive to periods of up to 10 days, with 
the same amplitude and for stars of the same mass. Thus, 
biasses could potentially have explained an deficit of slow 
rotators in Cep OB3b, but not an excess. 

Another possibility is contamination of our periodic 
sample. We have taken steps to remove foreground and back- 
ground contamination from our sample (see section l3~2Tl . but 
it remains possible that there is contamination from a sec- 
ond population of PMS stars towards Cepheus OB3b. This 
second population need not share the distance and age of 
the main population, so masses assigned to its members 
would be in error. In this way slowly rotating, high mass 
PMS stars might be confused for lower mass objects. Optical 
spectroscopy of the region could test this possibility. If sys- 
tematic error is not responsible for the observed difference 
between Cep OB3b and NGC 2362, it follows that differ- 
ent clusters at similar ages can show very different rotation 
distributions. This is not the first time such environmen- 
tal differences have been found; the young stars in IC348 
rotat e much slower than t hose in the similarly aged NGC 
2264 (|Littlefair et al. 20051 ) . 



Figure 13. Rotation period plotted against mass for the 
stars in Cepheus OB3b, as well as for other clusters with 
large collections of rotation dat within the litera ture. Ro- 
tation periods we re taken from iHerbst et al. (200211 (ONC) , 
lLamm et al. (200411 20051) (NGC 2264) and llrwin et al. (20081) 
(NGC 2362). See text for details of how masses were assigned. 



Cep OB3b. The difference is much more marked in other 
clusters; the same K-S test on the NGC 2362 sample gives a 
probability of less than 0.01 per cent that the high and low 
mass stars were drawn from the same distribution. Thus, 
whilst the stars in Cepheus show the same mass dependence 
of rotation seen elsewhere, the difference is less marked than 
is seen in other clusters. 



5.2 Rotation Period and Halpha emission 

In order to gauge the effect of accretion upon rotation 
rate we examine here the link betwe en rotation period 
and Hq emission. Ot her authors (e.g. iRebull et al. 20051 ; 
ICieza fc Baliber 20071 ) have examined the link between mid- 
infrared excess and rotation. Whilst an infrared excess is 
good evidence for a disc it is not evidence for ongoing ac- 
cretion, which is necessary for angular momentum regula- 
tion. Ha emission, on the other hand, is a signature of ongo- 
ing accretion, and mig ht also be expected to correlate with 
rotation rate. Indeed, lLamm et al. (20051 ) show that there 
is such a correlation in NGC 2264, in the sense that the 
stars showing strong Hq emission rotate more slowly, as a 
group, than those lacking Hq emission. A word of caution 
is needed: Hq emission strength is only an indicator of cur- 
rent accretion rate. Since the star-disc interaction, or a stel- 
lar wind, can only remove angular momentum on a finite 
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Figure 14. Period distributions for Cep OB3b, divided into sam- 
ples according to Ha strength. The stars are grouped in accreting, 
classical T-Tauri stars (CTTS) and non- accreting, Weak-lined T- 
Tauri stars (WTTS). See text for details of how stars were divided 
into these samples. 



timescale it would be more meaningful to compare rotation 
rate against the time- averaged accretion rate, taken over 
a comparable timescale. Additionally, it should always be 
born in mind that photometric monitoring can in principle 
be biassed against the detection of periods i n strongly ac- 
creting stars, although iLittlefair et al. (20051 ) showed that 
dense temporal sampling, as used in this paper, is effective 
in reducing such a bias. 

We divide our sample into accreting and non accret- 
ing objects in th e V — I, R — Ha diagram, following 
lLamm et al. (20041 ). Weak-line T-Tauri stars (WTTS) are 
assumed to be non-accreting, and are defined in this paper 
as objects with little or no Ha emission, as defined by 



(R - Ha) -(R- Ha) locus + A(R - Ha) < 0.0, 



(3) 



where (R — Ha)i ocus is defined by equation [T] Classical T- 
Tauri stars (CTTS) are objects with strong Ha emission, 
and are assumed to be accreting. In this paper, we define 
CTTS as those objects which satisfy 

(R - Ha) -{R- Ha) locu3 - A{R - Ha) > 0.1. (4) 

lLamm et al. (2005) show that equation 4 corresponds 
roughly to W\(H a ) > lOA, which is frequently used to select 
CTTS, based on the equivalent width of Ha. 

We find evidence that accretion is influencing rotation 
in Cepheus OB3b. The CTTS are rotating more slowly on 



average than the WTTS (see figure fTH This difference is sta- 
tistically significant; a ID K-S test gives a 5 per cent chance 
they were drawn from the same parent distribution. Thus, 
the data in Cepheus lends some support to the generally ac- 
cepted picture that accretion discs play a significant role in 
regulating the angular momentum of young stars. The result 
is only just statistically significant, however, and a statisti- 
cally significant difference between CTTS and WTTS is no 
longer present if the data is broken down into high- and 
low-mass sub-samples. This is perhaps not surprising. Our 
method of dividing stars into CTTS and WTTS is rather 
crude, relying as it does on narrowband photometry, rather 
than spectroscopy; it is not without sensitivity biasses which 
can be a function of mass. Also, we classify relatively few 
objects as CTTS, which reduces the statistical power of our 
sample. To test the influence of present-day accretion on 
rotation periods, good quality H-a spectroscopy of our pe- 
riodic objects is highly desirable. 



6 DISCUSSION 

Our results show that Cepheus OB3b is a rich cluster for 
studying the evolution of stellar angular momentum. The 
periodic data in Cep OB3b broadly confirms the general 
picture of angular momentum evolution in young stars; ac- 
creting objects are rotating more slowly than non-accretors, 
confirming the influence of discs on the angular momentum 
of young stars. As seen in other young clusters, the low mass 
stars are rotating more rapidly than the high mass stars, but 
this mass dependence is much weaker than that typically ob- 
served. Compared to the similarly aged cluster NGC 2362, 
Cep OB3b has very slowly rotating low-mass stars, and pos- 
sibly an excess of fast rotators amongst the high mass stars. 
There remains a possibility that this difference between the 
two clusters is a systematic error, perhaps related to con- 
tamination within our periodic database by a second PMS 
population. Optical spectroscopy of the region can test for 
this systematic effect, and can also improve our mass deter- 
minations, by determining individual extinction values for 
our periodic stars. 

If this difference between the two clusters is 
not systematic in origin, it confirms the findings of 
ILittlefair et al. (20051 ) , that clusters of similar ages can show 
very different rotational period distributions. Such environ- 
mental differences could possibly arise as a result of differ- 
ing angular momentum budgets between clusters, or they 
might reflect differences in disc lifetimes between star form- 
ing regions. These results should act as a cautionary tale 
for models which aim to reproduce the angular momentum 
evolution of stars (e.g. llrwin et al. 20071 ; iHerbst et al. 20021 ). 
These models are based on the assumption that the rotation 
period distributions of different clusters can be assembled 
into an evolutionary sequence, but this assumption is bro- 
ken if environmental differences between the clusters plays a 
larger role than evolutionary effects, due to magnetic wind 
braking, disc locking, or contraction towards the main se- 
quence. It is thus important to fit such models to datasets 
which include rotation periods from more than one cluster at 
each age, to gauge the impacts of environmental differences 
between clusters. 
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7 CONCLUSIONS 

We present a photometric study of / band variability to- 
wards the young association Cepheus OB3b. The study is 
sensitive to periodic variability on timescales of less than a 
day, to more than 20 days. The result is a database of 704 
periodic variables in the field of Cep OB3b. A random in- 
spection of 200 of these objects suggest that around 97 per 
cent of these periods are genuine. Colour cuts using V, I, R 
and narrowband Ha photometry are used to reject contam- 
inating objects, leaving 475 objects with measured rotation 
periods, which are very likely pre-main-sequence members 
of the Cep OB3b star forming region. 

We revise the distance and age to Cep OB3b, 
putting it on the con sistent age and distance ladder of 
iMavne fe Navlor ('20081 '). This yields a distance modulus of 
8.8±0.2 mags, corresponding to a distance of 580±60 pc, 
and an age of 4-5Myrs. For the purposes of this paper we 
therefore adopt an age of 4.5Myr. 

The rotation period distribution confirms the general 
picture of rotational evolution in young stars, exhibiting 
both the correlation between accretion and rotation ex- 
pected from disc locking, and the dependence of rotation 
upon mass that is seen in other star forming regions. How- 
ever, this mass dependence is much weaker than seen in 
other regions. Comparison to the similarly aged NGC 2362 
shows that the low-mass stars in Cep OB3b are rotating 
much more slowly. This points to a possible link between 
star forming environment and rotation properties. Such a 
link would call into question models of stellar angular mo- 
mentum evolution, which assume that associations can be 
assembled into an evolutionary sequence, thus ignoring en- 
vironmental effects. 
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